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Abstract: The last decade has seen considerable interest in the quasi-normal fre- 
quencies [QNFs] of black holes (and even wormholes) , both asymptotically flat and 
with cosmological horizons. There is wide agreement that the QNFs are often of 
the form Un = (offset) + in (gap), though some authors have encountered situations 
where this behaviour seems to fail. To get a better understanding of the general situ- 
ation we consider a semi-analytic model based on a piecewise Eckart (Poschl-Teller) 
potential, allowing for different heights and different rates of exponential falloff in 
the two asymptotic directions. This model is sufficiently general to capture and 
display key features of the black hole QNFs while simultaneously being analytically 
tractable, at least for asymptotically large imaginary parts of the QNFs. 
We shall derive an appropriate "quantization condition" for the asymptotic QNFs, 
and extract as much analytic information as possible. In particular, we shall explicitly 
verify that the (offset) + in (gap) behaviour is common but not universal, with this 
behaviour failing unless the ratio of rates of exponential falloff on the two sides of 
the potential is a rational number. (This is "common but not universal" in the sense 
that the rational numbers are dense in the reals.) We argue that this behaviour is 
likely to persist for black holes with cosmological horizons. 
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1. Introduction 



Black hole quasi-normal frequencies [QNFs] have been extensively studied over the 
past several decades, both for (potential) observational reasons, and for a number 
of highly technical theoretical reasons ^ |^, ||, |, ^. A key point is that the 
Regge-Wheeler equation, and the Zerelli equation, can (when expressed in terms of 
the tortoise coordinate) be written in the form 



drl 



[uj' -V{u)]i;. (1.1) 



The QNFs arising from this differential equation are in principle observable in astro- 
physical black hole oscillations and ring-down phenomena. At an abstract theoretical 
level it has been observed that the QNFs associated with this Regge-Wheeler equa- 
tion are very often of the form 

ojn = (offset) -|- in (gap). (1.2) 

Such behaviour was first noted when using the Eckart potential [|7| (Poschl-Teller 
potential [g, to estimate the lowest lying QNFs [|T2|, 0, g, but since then 
attention has focussed more on the higher lying QNFs as the mode number n — )■ 
oo [H, |15|, |l6l. For n — )■ oo the "gap" is simply related to the surface gravity of 



the black hole, at least in the asymptotically flat case where there is only a single 
horizon. The situation where both event and cosmological horizons are present is 
more subtle, and we shall seek a better semi-analytic understanding of this case by 
investigating a simplified model based on a piecewise Eckart potential (Poschl-Teller 
potential). This model has the virtues of being analytically tractable (at least for 
high lying QNFs as the mode number n — )■ oo) while being general enough to serve 
as a plausible model capturing the essence of the physics for a two-horizon black 
hole. 

We shall derive an appropriate "quantization condition" for the asymptotic 
QNFs, and from this condition extract as much analytic information as possible. 
In particular, we explicitly verify that the (offset) -|- in (gap) behaviour is common 
but not universal, with this behaviour failing unless the ratio of rates of exponential 
falloff on the two sides of the potential is a rational number. We argue that this 
behaviour is likely to persist for real asymptotically de Sitter black holes with both 
event and cosmological horizons. 

This (offset) (gap) behaviour is "common but not universal" in the sense that 
the rational numbers are dense in the reals. Furthermore, if one resorts to numerical 



^The nomenclature in this field is quite confusing, with the phrase "Poschl-Teller potential" 
being ascribed somewhat erratically to a number of different quantities. In this article we will 
always be referring to some variant of the sech^(-) potential, and will adopt the historically more 
accurate nomenclature of calling this (a special case of) the Eckart potential [Q. (Sec ^ for 
some general discussion of this historical point.) 



- 2 - 



techniques, then since the floating point numbers are rational approximations to real 
numbers, the (offset) + in (gap) behaviour is in fact universal within the context of 
floating point arithmetic, and all numerical experiments along these lines should be 
interpreted with this point in mind. 



2. Black hole QNFs 

For the specific case of a Schwarzschild black hole the tortoise coordinate is given by 



dr 



1 



2m 



r^(r) = r + 2m In 



r — 2m 
2m 



:2.ii 



and the Regge-Wheeler potential is 

. , / 2m\ 

V{n) = (i - — j 



£{£ + !) 2m(l - s^) 



(2.2) 



Here s is the spin of the particle and i is the angular momentum of the specific wave 
mode under consideration, with i > s. As r* — > — oo we have r — )■ 2m and 



V{r^) — )■ exp 



r* — 2m 
2m 



+ 1) + (1 - 
(2m)2 



Vq exp(2ftr*). 



(2.3) 



where k is the black hole surface gravity. This specific behaviour in terms of the 
surface gravity generalizes beyond the Schwarzschild black hole and for an arbitrary 
black hole in an asymptotically flat spacetime one has 



Vo- exp(— 2K|r*|), r* — )■ — oo; 
Vo+ (2m)2/r2, -> +oo. 



(2.4) 



This is often approximated (rather brutally) by an Eckart potential (Poschl-Teller 
potential) 

V{n) = J; (2.5) 
cosh [Kr^,) 

or more generally by using two adjustable parameters Vq and b 

Vo 



cosh^(r^,/6) 



(2.6) 



If one is interested in probing the lowest-lying QNFs one choses Vq and b to match 



the peak and peak curvature of the Regge-Wheeler potential 0, |T^, |Tj, [TS], |l| 



If instead one is interested in the higher QNFs one picks b = 1/k (and in this 
case chooses Vq in some more or less arbitrary manner) |T9|, |2^. Though 

the approximation is brutal, the Eckart potential has the great benefit of being 



3 



analytically soluble, and the QNFs are known exactly. (See for example |^|, or the 
discussion in 0. See also Taking 



\ - Vob'^ for 1/o&^ < 1/4, 



a 



(2.7) 



zJW-i for K)&'>l/4, 



the exact QNFs are 



iin + i ± a) 
^n,± = ^ 1 {n e {1, 2, 3, . . . }), 



(2.8) 



which is of the general form Un = (offset) + in (gap) with two families of QNFs 
described by the specific coefficients 



(offset) = 



(gap) = ^■ 



(2.9) 



Similar asymptotic behaviour for more realistic approximations to the Regge- Wheeler 
potential has been verified by a number of independent approaches. For Schwarzschild 
black holes, and more generally single-horizon black holes, it has been demonstrated 
that: 



WKB inspired approximations |]T^, |T3|, 



phase-amplitude methods [22, 23 



continued fraction approximations for the QNFs |]2^, p6| ; 
monodromy techniques |, |30 

Born approximations [0, [1^, |20| ; 



all agree on this general Un = (offset) + in (gap) form for the QNFs with the surface 
gravity n ^ 1/b governing the "gap" between the high-lying modes. For Reissner- 
Nordstrom, Kerr, Kerr-Newman, Sen [Reissner-Nordstrom-dilaton], black holes the 
situation is somewhat messier and sometimes still controversial, with results often 
being dimension dependent, and with results often (somewhat unexpectedly) de- 

31 El 



pending on the properties of the inner horizon |^ 



mm 



(See in particular for a recent review.) For example, even for the Reissner- 
Nordstrom black hole the Q/M ^ limit does not seem to reproduce the QNFs of 
the Schwarzschild black hole, and situations where -^/l — Q^/M^ is a rational number 
seem to be special. Though the techniques we adopt below are not directly applica- 
ble to the Reissner-Nordstrom black hole it is perhaps interesting to note that the 
special situations where ^1 — Q'^/M'^ is a rational number correspond to 
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the ratio of surface gravities at the inner and outer horizons, being the square of a 
rational number. 

If one turns to asymptotically de Sitter black holes (or more generally any two 
horizon system where the region of interest is bounded by the two horizons — an event 
horizon and a cosmological horizon) the situation is considerably more ambiguous — 
the Regge- Wheeler potential and the Zerelli potential have the asymptotic behaviour 

{Vo- exp(— 2K_|r*|), r* — )■ — oo; 
(2.10) 
^0+ exp(-2K+|r*|), r* -> +00; 

where the two surface gravities are now (in general) distinct. It is this situation that 
we will model using a piecewise Eckart potential (Poschl-Teller potential). A consid- 
erable amount of analytic information can be extracted from this model, information 
which in the concluding discussion we shall attempt to relate back to "realistic" black 
hole physics. 

3. Piecewise Eckart (Poschl-Teller) potential 
3.1 Potential 

The model we are interested in investigating is 

-V'"(x) + V{x) ipix) = 0, (3.1) 

with 

{Vo- sech^(a;/6_) for x < 0; 
(3.2) 
Vo+ sech^(a;/6+) for x > 0. 

We will allow a discontinuity in the potential at x = 0. The standard case that is 
usually dealt with is for 

Vo- = Vo+ = Vo; b- = b+ = b; V{x) = !° (3.3) 

cosh [x/b) 



A related model where Vo_ = Vo+ = Vq but 6+ 7^ 6_ has been explored by Suneeta [41 



but our current model is more general, and we will take the analysis much further. 
3.2 Wavefunction 

We start by imposing quasi-normal boundary conditions (outgoing radiation bound- 
ary conditions) 

V'+(a; ^ +00) ^ e"^^^; ^/'_(a; ^ -00) ^ e+''^^ (3.4) 
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On each half hne {x < 0, and x > 0) the exact wavefunction (see especially page 405 
of the article by Beyer is: 

1 



1 1 1 



where 



Vq\? for Vq\? < 1/4; 



a 



(3.5) 



(3.6) 



'Vob^ - \ for Voh'^ > 1/4. 
"All" we need to do is to appropriately match these wavefunctions at the origin. 

3.3 Junction condition 

The key step in matching these two wavefunctions at x = is to calculate the 
logarithmic derivative. Using the variable z = 1/(1 + e"*^^^/^*), note that a; = maps 
into z = 1/2. Then using the Leibnitz rule and the chain rule one has: 



1 din {2F1 (i + a±, i - a±, 1 + i6±cu, 2;)} 



^±(0) ' " ' 26± 
Invoking the differential identity 



dz 



(3.7) 



2=1/2. 



in appendix y, we see 



+iu 



2F1 (i + a±, I - «±, ib±u, z) 



2F1 (i + a±, \ - a±, 1 + ih±uo, z) 



(3. 



2=1/2. 



(In special situations where '?/'±(0) might accidentally equal zero one might need to 
perform a special case analysis. The generic situation is ^"^1(0) 7^ 0? ^"^^ will prove 
sufficient for almost everything we need to calculate.) Now using Bailey's theorem 
( |C.1| ) to evaluate the hypergeometric functions at 2; — )■ | we have the exact result 



V'±(o) 



_2_ r(^^^±^ + |)r( 
K Ti^^^±^ + i)r( 



3^ 

2 ' i' 

'a±^iu)b± I l^ 
2 1' 



(3.9) 



Now if u has a large positive imaginary part, then the Gamma function arguments 
above tend towards the negative real axis, a region where the Gamma function has 
many poles. This is computationally inconvenient, and to obtain a more tractable 
result it is extremely useful to use the reflection formula ( p.l|) of appendix to 
derive 

_ a±4^ _ 1^ sin(7r[ "^+r''^ + i)]) 



Yf a±+iu}b± _|_ 3^ 



' a±+iu)b± 



r(i 
rf- 



^^-l)sin(vr[^-''^ 
-|)sin(7r[ °^+^^"'^ 



a±+iu)b± I 3\ • / 

- 2 + 4) sin(vr 



■ + |]) 



r( 



.^±±^ + i) sin(7r[^^±:±^ + i] + f ! 

a±+iuib± I 3\ 



a-u + iuh-i- 1 



This leads to the exact resuh 

2 p/ ~a±-iuib± 



+ |)r 



■ a±-iujb± _|_ S-j 



&± r( 



X tan vr 



2 



1 \ a±—iujb± 
i)'- V 2 



I- iub-L. 1 
— +4 



tan TT 



-a± + iub± 1 
2 ^ 4 



(3.11) 



If w has a large positive imaginary part, then the Gamma function arguments above 
now tend towards the positive real axis, a region where the Gamma function is 
smoothly behaved — all potential poles in the logarithmic derivative have been 
isolated in the trigonometric functions. We can also use one of the trigonometric 
identities ( |A.1|) of appendix ^ to rewrite this as 



V4(o) 



2 Ti 



+ !)r( 



jwb-i- I 3^ 
2 4' 



cos(7ra;±) — cos(7r[ia;6± + 1/2]) 



+ i)r(^^±^ + i) 



which we can rewrite (still an exact result) as 
V^^(O) ^ r^-"±--f'± 



2 r( -"^r""^ + f; 



&± r( 



■p/a±-iwfe± _|_ 3^ 



jajfe± I 1 Np/ a±— , l^ 
2 4' ^ 2 4' 



cos(7ra-|-) + cos(7r[za;6± 



cos(7ra;±) + sin(i7ru;6±) 
cos(7rQ;-|-) — sin(i7ru;6±) 



1/2])' 
(3.12) 



(3.13) 



The exact junction condition we wish to apply at a; = is 

^;(o) ^-(0) 

^+(0) ^-(0)' 



(3.14) 



but the presence of the Gamma functions above makes this exact junction condition 
intractable. Fortunately, as long as we are primarily focussed on the highly damped 
QNFs (Im(c<;) — )■ oo) we can employ the Stirling approximation in the form ( |B.2|) 
indicated in appendix |B| to deduce 



p/ ±o.±-iixih± _|_ 3^ 



±a+ — %ijjh+ 1 

+ - X 

2 4 



Im(u;)6-i 



X 



1 + 



1 + 
1 

Im(u;6-| 



(3.15) 



This allows us to deduce an approximate junction condition for the asymptotic QNFs 
cos(7ra;+) + sin(z7rci;6+) cos(7rQ;_) + sin(27ra;6_) 



cos(7ra4.) — sin(i7ru;6_| 



cos(7ra_) — sinfivrwfe-' 



(3.16) 



which is accurate up to fractional corrections of order O (1/Im(ci;6±)). This is now 
an approximate "quantization condition" for calculating the QNFs, which is asymp- 
totically increasingly accurate for the highly-damped modes. 
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3.4 QNF condition 

The asymptotic QNF condition above can, by cross multiplication and the use of 
trigonometric identities, be rewritten in any one of the four equivalent forms: 



sin(— i7ra;6+) sin^—iTfcub-) = cos(7rQ;+) cos(7ra_); (3-17) 

sinh(7ra'6+) sinh(7ra;6_) = — cos(7rQ;+) cos(7rQ;_); (3.18) 

cos(— i7ra;[6+ — — cos{—i7iu[b^ + b_]) = 2 cos(7ra+) cos(7rQ;_); (3.19) 

cosh(7ra;[6+ — 6_]) — cosh(7ra;[6+ + 6_]) = 2 cos(7rQ!+) cos(7rQ;_). (3.20) 



Which particular form one chooses to use is a matter of taste that depends on exactly 
what one is trying to establish. It is sometimes useful to split the asymptotic QNF 
condition into real and imaginary parts. To do so note 

cos{A + iB) — cos{A) cos{iB) — sm{A) sm{iB) 

= cos(A) cosh(5) - i sm{A) sinh(5), (3.21) 

so that 

cos{—i7iu![b^ — 6_]) = cos(lm{u)n[b+ — 6_]) cosh(Re(w)7r[6+ — 6_]) 

-isin(Im(6j)7r|6+ - 6_|) sinh(Re(6j)7r|6+ - 6_|). (3.22) 

Therefore the asymptotic QNF condition imphes both 

cos(Im(a;)7r[6+ — 6_]) cosh(Re(co')7r[6+ — 6_]) 

= cos(Im(a;)7r[6_|_ + 6_]) cosh(Re(a;)7r[6+ + 6_]) + 2 cos(7rQ;+) cos(7rQ;_), (3.23) 

and 

sin(Im(a;)7r|6+ — 6_|) sinh(Re(a;)7r|6+ — 6_|) 

= sin(Im(a;)7r[6+ + 6_]) sinh(Re(a;)7r[6+ + 6_]). (3.24) 

We shall now seek to apply this QNF condition, in its many equivalent forms, to 
extract as much information as possible regarding the distribution of the QNFs. 

4. Some general observations 

We shall start with some general observations regarding the QNFs. 

1. Note that the a± are either pure real or pure imaginary. 

2. Consequently cos(7rQ;+) cos(7rQ;_) is always pure real e [— l,+oo). 
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3. If cos(7ra+) cos(7ra_) > 0, then there are no pure real QNFs. 



Proof: Consider equation ( p.l8| ) and note that under this condition the LHS is 



positive while the RHS is negative. 
4. If cos(7ra;_|_) cos(7ra_) < 0, then there is a pure real QNF. 



Proof: Consider equation ( p.l8 ) and note that under this condition the RHS 



is positive. The LHS is positive and by continuity there will be a real root 
CO G (0, oo). 

5. If cos(7ra+) cos(7ra_) > 1, then there are no pure imaginary QNFs. 



Proof: Consider equation (|3.17| ) and note that under this condition the LHS 
< 1 while the RHS > 1. 

6. There are infinitely many pure imaginary solutions to these asymptotic QNF 
conditions provided cos(7ra+) cos(7ra_) < Q{b^,b_) < 1; that is, whenever 
cos(7ra+) cos(7ra_) is "sufficiently far" below 1. 

Proof: Define 



Then by inspection equation (|3.19| ) will have an infinite number of pure imag- 
inary solutions as long as 

cos(7rQ;+) cos(7ra_) < Q{b+,b_). (4.2) 

7. For any purely imaginary uj there will be some choice of b±, a± that makes 
this a solution of the asymptotic QNF condition. 



Proof: Consider the specific case 

^ ^ i{a . ^ ^^^^^ 



and note this satisfies the QNF condition but enforces only two constraints 
among the four unknowns b±, Vo±. 

5. Some approximate results 

A number of approximate results can be extracted by looking at special regions of 
parameter space. 
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5.1 Ceise 6_ ^ b+ 

Suppose 6_ ^ b^. Then the quantity —i{b^ — b-.)uj is slowly varying over the range 
where —i{b^ + b_)uj changes by 27r. Let be any solution of the approximate QNF 
condition, and define 6* = (6+ + b-)/2. Then for nearby frequencies we are trying to 
(approximately) solve 

cos(— ia;*7r[6+ — — cos(— ia'7r[6+ + b-]) — 2 cos(7rQ!+) cos(7rQ;_), (5.1) 

that is 

cos(— ia;*7r[6+ — 6_]) — cos(— i<x;27r6*) = 2cos(7rQ;+) cos(7rQ;_), (5.2) 

and the solutions of this are approximately 

in , „ , , 6+ + 6- . 
~ a;* + — valid for \n\ <S ■- ; — -. (5.3) 

Thus approximate result will subsequently be incorporated into a more general per- 
turbative result to be discussed below. 

5.2 Case 6_ < 6+ 

Now suppose b- <^ 6+. Then the quantity —ib-ou is slowly varying over the range 
where —ib+cu changes by 2n. Let a;* be any solution of the approximate QNF condi- 
tion, then for nearby frequencies we are trying to (approximately) solve 

sm{—iuj7rb+) sin(— ia;*7r6_) = cos(7rQ;+) cos(7rQ;_), (5.4) 

and the solutions of this are approximately 

cun ^ cj* + — valid for \n\ <C — . (5.5) 

5.3 Case «_ fa 1/2 

This corresponds to 

Vo-bt ^ 0, (5.6) 
in which case the QNF condition becomes 

sin(— ia;7r6+) sin(— ia;7r6_) ~ 0. (5.7) 

Therefore one obtains either (the physically relevant condition) 

in 

—iojb+ — n =^ oj — — , (5.8) 

b+ 

or (the physically uninteresting situation) 

—iuib- = n. (5.9) 

Note: If you go to the limit a_ = 1/2 by setting Vo_ = 0, one sees on physical grounds 
that 6_ is irrelevant, so it cannot contribute to the physical QNF. Alternatively if 
you hold Vq- ^ but drive 6_ — > 0, then these QNF's are driven to infinity — and 
so decouple from the physics. Either way, the only physically interesting QNFs are 
uj — in/b+. We explore these limits more fully below. 
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6. Some special cases 

A number of special cases can now be analyzed in detail to give us an overall feel for 
the general situation. 

6.1 Case «_ = = «+ 
This corresponds to 

Vo-b'_ = ^ = Vo+bl, (6.1) 
in which case the QNF condition becomes 

sin(ia;7r6+) sin(ic(;7r6_) = 1. (6.2) 

Let us look for pure imaginary QNFs. (We do not claim that these are the only 
QNFs.) This implies that we must simultaneously satisfy both 



or both 
That is both 
or both 

Therefore either 



sm{—iuj7rb+) — sm{—iuj7rb-) — 1, (6.3) 

sm{—iu!7ib+) = sm{—iu!7ib_) = —1. (6.4) 

-ia;6+ = 2n+ + -; —iub- = 2n_ + -, (6.5) 

1 1 

-icub^ = 2n+ — -; —iub_ — 2n_ — -. (6.6) 



b+ ^ 2n+ + \ ^ An+ + 1 b^ ^ 2n+ - \ ^ 4n+ - 1 

6_ 2n_ + 1 4n_ + 1 ' 6_ 2n_ - | 4n_ - 1 ' ^ ' ^ 

In either case we need to be rational, so that 6+ = and 6_ = This 

special case is thus evidence that there is something very special about the situation 
where b+/b- is rational, more on this point below. Then either 

—ILO = ; -] —lU = ; -] (6.8) 



or 



Now write 



2n+-\ 2n_ - ^ ^ ^ 



n+ = m+ + n_ = m_ + np_; (6.10) 

with m+ < p+ and m_ < p_. (While n e {0, 1, 2, 3, . . . }.) Then in the first case 

r2m+ + | n1 (2m- + l nl m fell) 
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while in the second case 

. ( 2m+ - I n 1 . r 2m_ - | n 1 in 
\ p+K K j \ P-b* K j * K 

This is a precursor of the much more general result that we shall ultimately obtain 
for generic rational 

6.2 Case Vq^ = 

We can best analyze this situation by working directly with the exact wavefunction. 
If Vo- = then a_ = 1/2 and 

V'-(0) = 1; V'-(0) = +ia;; ^^^^ ^ +'^- (6-13) 
The exact QNF boundary condition is then 

= -r- ] !■ , f . (6.14) 

But this we can rewrite as 

^ ^ 2^ r(^^^^ + |)r(^^±^ + I) ^ cos(7ra+)-sin(-ma;&+) ^ 
6+ p/^a+-ic£6+ _^ iNp/a+-it£6+ _^ 1) ^ cos(7rQ;+) + sin(-i7ra;6+) ' ^ ' 

This certainly has pure imaginary roots. If we write ui — then asymptotically 
(jo; I — >■ oo) this becomes 

cos(7rQ;+) — sin(7r|c(;|6+) 



cos(7rQ;+) + sin(7r|a'|6+) ' 



(6.16) 



implying 



sin(7r|a;|6+) = 0; =^ 7r|a;|6+ = wk] =^ u) = —. (6.17) 

This agrees with our previous calculation for «_ Ri 1/2. 
6.3 Case b- ^ 

This is best dealt with by using a Taylor expansion to show that 

"^-^^^ +iuj + Vo-b. + 0{b''_). (6.18) 



V^_(0) 
That is 

^' (0) . . 

lim = 6.19 

The analysis then follows that for the case Vo^ — above, and furthermore agrees 
with our previous calculation for «_ ?s 1/2. 
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6.4 Ceise 6_ ^ oo 

This is best dealt with by using the Stirhng approximation together with a Taylor 
expansion to show that 

That is 



iy/iu^-Vo- + 0{l/bl). (6.20) 



hm ^ = t^u;^-Vo.. (6.21) 



The exact QNF boundary condition is then 



V p(a+:^ ^ l)r( -"+y"''+ + i) ■ 



(6.22) 



But this we can rewrite as 



_ ^ 2^ r(^i^±^ + |)r(^^±^ + I) ^ cos(7ra+) - sin(-z7ra;6 



p^-a+-ia^ _^ l-jp^a+-ia;fe+ _^ 1^ COs(7rQ;+) + Sm{-iTTU;b+) ' 

^ ^ (6.23) 

If we write cu = i\u!\ then asymptotically, (|a;| — > oo, with Vq- held fixed, implying 
that Vo- efTectively decouples from the calculation), this becomes 

^ ^ cos(7ra+) -sin(7r|a;|6+) 
cos(7rQ;+) + sin(7r|(A;|6+) ' 

implying 

in 

sin (tt leu 1 6+) = 0; =^ 7r|c(;|6+ = nir; =^ cu — —. (6.25) 

b+ 

The importance of this observation is that it indicates that for "one sided" potentials 
it is only the side for which the potential has exponential falloff that contributes to 
the "gap". 



7. Rational and irrational ratios for the falloff 

We have already seen above good reason to suspect that rational ratios might 
be special. Let us now explore this case in more detail. 

7.1 Explicit examples 

• If 6+ = 6_ = 6*, but we do not necessarily demand = then the asymp- 
totic QNFs are exactly calculable and are given by 

i cos~-^ {1 — 2 cos(7rQ;+) cos(7ra_)| in 
= 2.6. + T- (^-1) 
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Proof: The asymptotic QNF condition reduces to 



1 — cos(— i27ra;6*) = 2 cos(7rQ;_|_) cos(7rci;_), (7.2) 

whence 

cos{—i2'KUjb^) = 1 — 2 cos(7rQ;+) cos(7rQ;_). (7.3) 
This is easily solved to yield 

—i2Tiu}nK = cos~^ {1 ~ 2 cos(7rQ;_|_) cos(7rQ;_)} + 2mT, (7.4) 

whence 

icos~^ {1 — 2 cos(7rQ;+) cos(7rQ;_)} ^ in 
2-nh^ 6* 

Comment: These QNFs are pure imaginary for cos(7ra+) cos(7rQ;_) < 1, and 
off-axis complex for cos(7ra+) cos(7ra_) > 1. We can always, for convenience, 
choose to define Re(cos~^(x)) e [0, 27r); then cos~^(-) is double valued. 

9 e [0, tt] and cos^ = x implies cos~^(a;) = {^, tt — 9}, (7.6) 

9 e [tt, 27r) and cos^ — x implies cos~^(a;) = {9, Sn — 9}. (7.7) 

With this notation 

in icos~-^{l — 2 cos(7rQ;+) cos(7rQ;_)} 
ujn — ujo + 7—; ujo — r-^^ ; (7.8) 

with < Im(cjo) < V^* ^ ^ {0,1,2,3,...}. Note that because of the 
double valued nature of cos~^(-) there are actually two branches of QNFs hiding 
in this notation — which we will need if we wish to regain the known standard 
result when we specialize to q;_ = «+. (We shall subsequently generalize 
this specific result, but it is explicit enough and compact enough to make it 
worthwhile presenting it in full. Furthermore we shall need this as input to our 
perturbative analysis.) 

If 6+ = 36- , that is 6+ = and &_ = ^b^, but we do not necessarily demand 
a+ — a-, then the asymptotic QNFs are calculable and are given by 



i 1 / 1 ± a/Q — 16 cos(7rQ;+) cos(7rQ;_)) \ 2m 



Proof: To see this note that in this situation — 6_| =6* = (6+ + h-)/2. 
Therefore the QNF condition reduces to 

cos(— i7ra'6*) — cos(— i27ra;6*) = 2 cos(7rQ;+) cos(7rQ;_), (7-10) 
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implying 

cos(— i7ra;6*) — 2 cos^ {—inub^) + 1 — 2 cos(7rQ;+) cos(7rQ;_). (7-11) 
That is 

2 cos^(— i7ra;6*) — cos{—iTTUjb^) — 1 + 2 cos(7ra+) cos(7rQ;_) = 0, (7-12) 
whence 

1 ± a/I + 8(1 — 2 cos(7rQ;+) cosina-)) 

cos{-i7rujb^) = ^ — -, (7.13) 

so that 

1 ± \/9 — 16 cos(7rQ;+) cos(7rQ;_)) 

cosi-iTTCub^) = ^ ^ — — ^ -, (7.14) 

implying 

1 / 1 dz a/9 - 16 cos(7rQ;+) cos(7rQ;_)) , ^ . 

-inunb* = cos"^ | ^ j — ±^ ^ ^ ) + n2n. (7.15) 

Finally 



i _ | / 1 ± -4/9 — 16 cos(7rQ;+) cos(7rQ;_)) \ 2m 

= rt: 4 j + X- 

Comment: This gives us another specific example of asymptotic off-axis com- 
plex QNF's — now with bj^^b^. Note that because of the ± and the double- 
vahicd nature of cos~^(-) there are actually 4 branches of QNFs hiding in this 
notation. 

• This particular trick can certainly be extended to the cubic and quartic poly- 
nomials, for which general solutions exist. 



The quadratic corresponds to 

b^-b. 



b+ + b. 

The cubic corresponds to 

b^-b. 



b+ + b_ 

- The quartic corresponds to 

6+ - &- . , 5 



3 6_. (7.17) 



3; 6+ = 2 6_. (7.18) 



= 4; b+ = - 6_. (7.19) 



+ ' "3 



This very strongly suggests there is something special about arbitrary rational 
values of In fact the current results will very shortly be subsumed into 

much more general results for rational values of 
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7.2 Rational ratios for the fallofF 

Suppose b+/b- is rational, that is 



b- p- ' 



(7.20) 



and suppose we now define 6* by 

b^ ^ p^b^; 6* = hcf(6+, (7-21) 

then the asymptotic QNF condition is given by 

sin(— ia;7rp+6*) sm{—iojnp_b^) = cos(7rQ;+) cos(7rQ;_). (7.22) 

If cj* is any specific solution of this equation, then (setting s = 1 for p+p- odd, and 
s = 2 for p+p- even) 

ins / 1 1 \ 

ujn = <^* + -r- = + ins 1cm I — , — I (7.23) 



6* ^ 

will also be a solution. But are these the only solutions? Most definitely not. For 
instance, consider (for rational &+/&-) the set of all QNFs for which 

Im(a;) < ^, (7.24) 

and label them as 

uJo,a a e {1,2, 3... TV}. (7.25) 
Then the set of all QNFs decomposes into a set of famihes 

ZTIS 

0On,a^uJo,a+-^; a G {1, 2, 3 . . . TV}; n G {0, 1, 2, 3 . . . }; (7.26) 

where is yet to be determined. But for rational b^/b^ we can rewrite the QNF 
condition as 

cos(— ia;7r6*|p+ — p_|) — cos(— ia;7r6*[p+ +p_]) = 2 cos (7rQ;+) cos (ttq;-) . (7.27) 
Now define z = exp(a;7r6*), then the QNF condition can be rewritten as 

z\p+-P-\ + ^-|P+-P-I _ ^b++P-] _ z-\P++P-] ^ 4 cos {7ra+) cos (ttq;.) , (7.28) 

or equivalently 

z-'\p++p-]-z^p+-p-\+\p++p-]+4^cos{na+) cos(7rQ;_) ^+b++P-] _^-|p+-p-l+b++P-]+i ^ q, 

(7.29) 



-16- 



that is 

^2[p++p_] _ ^2p^,. ^ 4 ^+b++P-] _ ^2p^i„ + 1 = 0. (7.30) 

This is a polynomial of degree = 2(p++p_), so it has exactly N roots Za (occurring 
in complex conjugate pairs, and whenever is odd, in symmetric pairs ±Za). 

Then the QNFs are, with the imaginary part of the logarithm lying in [0, 2tt), 

ini Z ) Ills 

^n,a = ^j^ + -r- a e {1,2, 3... s{p++p_)} n e {0,1,2,3, .. .}. (7.31) 

So for rational b^/b_ with b^/b_ = p+/p- we have exactly s(p+ + P-) equi-spaced 
families of QNFs all with with gap and with (typically distinct) offsets \n{za) / (7r6*). 
That is: Arbitrary rational ratios ofb^/b- automatically imply the Un = (offset) + 
in (gap) behaviour. 

Comparison with some monodromy results: Consider the specific situation 
presented in that is, Schwarzschild-de Sitter (Kottler) spacetime with surface 
gravity kn at the black hole horizon and kc at the cosmological horizon. In that 
article monodromy techniques are used to obtain an approximate QNF condition (for 
scalar fields) of the form 

cosh ^ + 3 cosh — + — = 0. (7.32) 

\kH kc J \kH kc J 

Qualitatively similar QNF conditions hold for higher spin. 

Now if kn/kc = nn/nc is rational then we can (strategically) choose integers 
nc and uh so that ku/kc = nc/nH, and so write ku = k^jriH while kc = k^^jnc- 
Then this approximate QNF condition becomes 

eosh [ ^ - + 3 cosh [ + ^^'j = 0. (7.33) 

Defining z = exp(7ra;/A;*) this becomes the (Laurent) polynomial condition 

which we can rearrange as the (ordinary) polynomial 

3^2(nH+nc) ^ ^2nH ^ ^2nc +3 = 9. (7.35) 

This is an ordinary polynomial with degree = 2{nH + nc) so it has = 2{nH+nc) 
roots Za (typically distinct, apart from the fact that they occur in complex conjugate 
and symmetric pairs i^a)- Thus in a manner qualitatively similar to our direct 
semi-analytic results above, the monodromy based approximate QNFs will fall into 
families of the form 

In ( z ) tn 

uJn,a = ^7^ + Tr a & {1,2,3 ... {uh + Uc)} n G {0, 1, 2, 3, . . . }. (7.36) 
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In each individual family the real part of the QNF does not oscillate. As one moves 
from one family to the next the real part of the QNF will change. If the degree of the 
polynomial N = 2{nH + nc) is sufficiently high then the discrete "base frequencies" 
hi{za) / Trk^ will visually approximate a single cycle of an oscillating curve. In figure 3 
of ref several thousand famihes of equi-spaced QNFs combine to yield a discrete 
collection of points that visually appears to be a repetitive oscillating curve. That 
is, the present semi-analytic techniques, where they overlap, are compatible with 
expectations based on monodromy techniques. 

7.3 Irrational ratios for the falloff 

Now suppose is irrational, that is 

6, = hcf(6+,6_) = 0. (7.37) 

Then all of the "families" considered above only have one element 

cuca a G {l,2,3...cx)}. (7.38) 

That is, there will be no "pattern" in the QNFs, and they will not be regularly 
spaced. (Conversely, if there is a "pattern" then is rational.) Stated more 

formally, it is possible to derive a theorem as below. 

Theorem: Suppose we have at least one family of equi-spaced QNFs such that 

ujn = + inK, (7.39) 

then is rational. 



Proof: If we have a family of QNFs of the form given in equation (|7.39|) then we 
know that \/n > 

cos{—iuo'K\b^ — b-\ + nKn\b+ — 6_|) — cos(— 'i(X'o7r[6+ + 6_] + n_ft'7r|6+ + 6_|) 

= cos{—iijJo7i\b+ — 6_|) — cos(— 'i(X'o7r[6+ + 6_]). (7.40) 

Let us write this in the form Vra > 

cos{A + nJ) - cos(5 + nL) = cos(A) - cos(5), (7.41) 

and realize that this also implies 

cos(A +[n + l]J)- cos(5 + [n + 1]L) = cos(A) - cos(5), (7.42) 

and 

cos(A +[n + 2]J)- cos(B + [n + 2]L) = cos(A) - cos(5). (7.43) 
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Now appeal to the trigonometric identity (based on equation ( |A.3|) ) 

cos(A +[n + 2]J) + cos{A + nJ) = 2 cos( J) cos(A + [n + 1] J), (7.44) 
to deduce 

cos(J) cos(A +[n + 1] J) - cos(L) cos(5 +[n + l]L) = cos(A) - cos(5). (7.45) 
That is, > we have both 

cos(A + [n + 1] J) - cos(5 + [n + 1]L) = cos(A) - cos(5), (7.46) 

and 

cos(J) cos(v4 + [n + 1] J) - cos(L) cos(5 + [ra + 1]L) = cos(v4) - cos(5). (7.47) 

The first of these equations asserts that all the points 

cos(A +[n + l]J), cos(5 + + 1]L)) (7.48) 

lie on the straight line of slope 1 that passes through the point (O,cosi? — cos A). 
The second of these equations asserts that all the points 

cos(A +[n + l]J), cos(5 + [n + 1]L)) (7.49) 

also lie on the straight line of slope cos(J)/ cos(L) that passes through the point 
(0, [cos-B — cosA]/ cosL). We then argue as follows: 

• If cos J 7^ cos L then these two lines are not parallel and so meet only at a 
single point, let's call it (cosA*, cos-B*), whence we deduce 

cos{A+[n + l]J) = cosA^; cos(5 + [n + 1]L) = cos S*. (7.50) 

But then both J and L must be multiples of 27r, and so cos J = 1 = cosL 
contrary to hypothesis. 

• If cos J = cos L 1 then we have both 

cos(A +[n + 1]J) - cos(fi +[n + 1]L) = cos(A) - cos(5), (7.51) 

and 

cos(J) [cos(A +[n + 1] J) - cos(fi + [n + 1]L)] = cos(A) - cos(5). (7.52) 

but these are two parallel lines, both of slope 1, that never intersect unless 
cos(J) = 1. Thus cos J = 1 = cosL contrary to hypothesis. 
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Wc therefore conclude that both J and L must be muhiples of 2ti so that 
cos J = 1 = cos L (in which case the QNF condition is certainly satisfied) . 



But now 

\h+-h-\ J 



- e Q, (7.53) 



b+ + b- L 
and therefore 

^ e g. (7.54) 

That is: Rational ratios ofb+/b- are implied by the cun — (offset) +in (gap) behaviour. 

8. Systematic first-order perturbation theory 

Sometimes it is worthwhile to adopt a perturbative approach and to estimate shifts 
in the QNFs from some idealized pattern. Define 

b^h±^. A=|6+-6_|; (8.1) 
and rewrite the asymptotic QNF condition as 

cos(— i7ra;A) — cos(— i27ra;6) = 2 cos(7rQ;+) cos(7rQ;_), (8.2) 

where we are implicitly holding a± fixed. When A = we have previously seen that 
the QNF are exphcitly calculable with 

i cos"-^ {1 — 2 cos(7rQ;+) cos(7rQ;_)| in ,^ „, 

= 2nb +T <«-^) 

Can we now obtain an approximate formula for the the QNF's when A 7^ 0? It is a 
good strategy to define the dimensionless parameter e by 

A = 2 e 6, (8.4) 

and to set 

uj — uj + Sou; Sou — 0{e); (8.5) 
so that the asymptotic QNF condition becomes 

cos(— i27r[c<; + Suj]eb) — cos(— i27r[a; + 6uj]b) = 2 cos(7rQ;_|_) cos(7rQ;_). (8.6) 

Then to first order in e 

cos(— i27ra)e6) — cos(— i27r[a; + 6u!]b) — 2 cos(7ra+) cos(7ra_), (8.7) 

where implicitly this approximation requires e|(5a;|6 -C 1. Subject to this condition 
we have 

cos(— i27r[a; + Su!]b) — cos(— i27ra;e6) — 2 cos(7rQ;+) cos(7rQ;_), (8.8) 
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whence 

—i2TT[uJn + 5un]b = cos^^ {cos(— i27rci;„e6) — 2 cos(7rQ;+) cos(7rQ;_)} + 2™, (8.9) 
so that 

.cos^""^ {cos(— z27r(i;„e6) — 2 cos(7rQ;+) cos(7ra_)} in /o inN 

UJn + dUn = i h — . (8.10) 

Zno 

But we know that the unperturbed QNFs satisfy Un = uq + in/b, so we can also 
write this as 

.cos^^ {cos(— 'i27rw„e6) — 2 cos(7ra+) cos(7ra_)} ^ /o tt\ 

dUn = I — Wo- (o-ll) 

Itto 

Using the definition of Con this can now be cast in the form 

.cos~^ {cos(— i27ra;„e6) + cos(— i27rw„6) — 1} „ /o i o\ 

OUJn = 1 —J ^0, (o-l^) 

/TTO 

or the shghtly more suggestive 

.cos~^ {cos(-i27rw„6) + cos(-i27rc2;„e6) - 1} ^ /o i o\ 

OUn = % —-, Wo, {6.16} 

Z-KO 

which can even be simphfied to 

. cos~^ {cos(— 'i27rwo&) + cos(— i27rci;„e6) — 1} „ /o i /i\ 

OUn = i ^^—^ Wo. (o-14) 

Note that this manifestly has the correct limit as e — )■ 0. Note that we have not 
asserted or required that (W„e6 ^ 1, in fact when n ^ 1/e this is typically not true. 
(Consequently cos(— i27rc2;„e6) is relatively unconstrained.) Note furthermore that 
Im(5c<;„) < 1/6. 



9. Discussion 



The key lesson to be learned from our semi-analytic model for the QNFs is that 
the Un = (offset) + in (gap) behaviour is common hut not universal. Specifically, in 
our semi-analytic model the key point is whether or not the ratio is a rational 

number. Let us also note here that monodromy techniques (see in particular [^j) quite 
often also lead to the qualitative result that rational ratios of surface gravities are 
closely related to the occurence of w„ = (offset) -|-m (gap) behaviour. This behaviour 
is thus "common but not universal" in the sense that mathematically the rational 
numbers are dense in the reals, while when applying numerical techniques fioating 
point numbers are a subset of the rational numbers. So the w^ = (offset) + in (gap) 
behaviour is actually universal within the context of fioating point arithmetic. 
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Will this behaviour extend to more "realistic" astrophysical or asymptotically 
de Sitter black holes? Consider a "wavepacket" , built up out of highly damped modes, 
that is centered near the peak of the Regge- Wheeler (Zerelli) potential. While the 
initial short-time behaviour of the wavepacket is likely to be sensitive to the details of 
the Regge-Wheeler (Zerelli) potential, such a wavepacket will quickly damp out and 
spread out towards both r* — >■ — oo and r* — > +00, so that the wavepacket will pen- 
etrate regions where our piecewise Eckart model potential, (governed by the surface 
gravities at the event and cosmological horizons), should be a good approximation 
to the true potential. We should therefore expect the results of our semi-analytic 
model to be qualitatively (but not necessarily quantitatively) accurate for estimating 
the asymptotic QNFs of "realistic" asymptotically dc Sitter black holes. Because of 
the way the asymptotic QNF condition was derived, we do not expect out model to 
give good results for low-lying QNFs. 

Overall, one of the nice features of this semi-analytic model is that a quite sur- 
prising amount of semi-analytic information can be extracted, in terms of general 
qualitative results, approximate results, perturbative results, and reasonably explicit 
computations. We suspect that it might be possible to generalize the model potential 
even further — the "art" would lie in picking a piecewise potential that is still ana- 
lytically solvable (at least for the highly damped modes) but which might be closer in 
spirit to the Regge-Wheeler (Zerelli) potential that is the key physical motivation for 
the current article. (Of course if we temporarily forget the black hole motivation, it 
may already be of some mathematical and physical interest that we have a nontrivial 
extension of the Eckart potential that is asymptotically exactly solvable — one could 
in principle loop back to Eckart 's original article and start asking questions about 
tunelling probabilities for electrons encountering such piecewise Eckart barriers.) 

A. Trigonometric identities 

In the body of the article we will need to use some slightly unusual trigonometric 
identities. They can be derived from standard ones without too much difficulty but 
are sufficiently unusual to be worth mentioning explicitly: 



tan A tan B — 



cos(A - B) - cos(A + B) 



(A.l) 



cos(^ - B)+ cos{A + B) 




(A.2) 



and 



cos( A 2B) cos A = 2 cos B cos{A + B) . 



(A.3) 
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B. Gamma function results 



The key Gamma function identity we need is 

r(^) r(i - ^) = (B.i) 

We also need the foUowing asymptotic estimate based on the Stirhng approximation 

r(^ + l)_ 



z 



Re(;2) ^ oo. (B.2) 



V{z) 

C. Hypergeometric function identities 

The key hypergeometric function identities we need are Bailey's theorem 

which is easily found in many standard references, and the particular differential 
identity 

d {2^1(0, 6, c,^)} c-1 T? ( u w {no\ 

= [ 2F1 (a, b,c-l,z)- 2F1 (a, 6, c, z)\ , (C.2) 

az z 

which is not found in any of the standard references (but is easy enough to verify 
once it has been presented). 
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